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ABSTRACT 

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have assembled a complete sample of 
20 iGselected galaxies at 2.0 < z < 2.7 with high quality near-infrared spectra. As described in a 
previous paper, 9 of these 20 galaxies have strongly suppressed star formation and no detected emission 
lines. The present paper concerns the 11 galaxies with detected Ha emission, and studies the origin 
of the line emission using the GNIRS spectra and follow-up observations with SINFONI on the VLT. 
Based on their [Nil] /Ha ratios, the spatial extent of the line emission and several other diagnostics, 
we infer that four of the eleven emission-line galaxies host narrow line active galactic nuclei (AGNs). 
The AGN host galaxies have stellar populations ranging from evolved to star-forming. Combining 
our sample with a UV-selected galaxy sample at the same redshift that spans a broader range in 
stellar mass, we find that black-hole accretion is more effective at the high-mass end of the galaxy 
distribution (~ 2.9 x 1Q 11 Mq) at z ~ 2.3. Furthermore, by comparing our results with SDSS data, we 
show that the AGN activity in massive galaxies has decreased significantly between z ~ 2.3 and z ~ 0. 
AGNs with similar normalized accretion rates as those detected in our if-selected galaxies reside in 
less massive galaxies (~ 4.0 x 10 10 M©) at low redshift. This is direct evidence for downsizing of 
AGN host galaxies. Finally, we speculate that the typical stellar mass-scale of the actively accreting 
AGN host galaxies, both at low and at high redshift, might be similar to the mass-scale at which 
star-forming galaxies seem to transform into red, passive systems. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift 
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1. INTRODUCTION 

The galaxy population today can very broadly be di- 
vided into a population of star-forming disk galaxies and 
a population of passive early-type galaxies. Although 
this has been known for a long time, the physics behind 
the dichotomy in galaxy properties are still poorly un- 
derstood. In particular, the model of the formation of 
massive, red elliptical galaxies has frequently been mod- 
ified, driven by observations. The recent popular hierar- 
chical galaxy formation models produce galaxies of the 
required mass, but can not match their quiescent st el- 



Electronic address: mariska@strw.leidenuniv.nl 

1 Based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile (076.A-0464 and 076.A-0718) 

2 Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership. 

3 Leiden Observatory, Leiden University, PO Box 9513, 2300 RA 
Leiden, The Netherlands 

4 Department of Astronomy, Yale University, New Haven, CT 
06520 

5 Yale Center for Astronomy and Astrophysics, Yale University, 
New Haven, CT 06520 

6 UCO/Lick Observatory, University of California, Santa Cruz, 
CA 95064 

7 Max-Planck-Institut fur extraterrestrische Physik, Giessen- 
bachstrasse, Postfach 1312, D-85748 Garching, Germany 

8 NSF Astronomy and Astrophysics Postdoctoral Fellow 

9 Carnegie Fellow, Carnegie Observatories, 813 Santa Barbara 
Street, Pasadena, CA 91101 

10 Departamento de Astronomi'a, Universidad de Chile, Casilla 
36-D, Santiago, Chile 

Goldberg Fellow, National Optical Astronomy Observatory, 
950 North Cherry Avenue, Tucson, AZ 85719 



lar populations, unless a mechanism - such as feedback 
from active galactic nuclei (AGNs) - is invoked to stop 
star formation at early times (e.g., |G ranato et al J 12004 ; 
ICroton et alll2006t iBower et al.ll2006t iKang et alJ[2006t ). 

Understanding the role of AGNs and in particular their 
feedback processes in the star formation history (SFH) of 
a galaxy is one of today's major challenges. The tight re- 
lation b etween the black hole mass and bulge velocity dis- 
persion (jFerrarese fc MerrittfeOOOHGebhardt et al.ll2000h 
may imply that black hole accretion is directly related to 
the formation of its host galaxy. Moreover, stellar popu- 
lations of the hos t galaxies seem to relate to the strength 
of the AGN (e.g. JKauffmann et al.l[2003bl ). However, the 
effects of AGN feedback processes are still poorly under- 
stood. The best example of AGN feedback "at work" is 
the brightest cluster galaxy Perseus A, w hich is inject- 
ing energy in the intracluster medium (e.g.. lFabian et ail 
120031 l2006f ) . There are several other recent attempts to 
constrain the truncation mechanism, which use the large 
stat istical data-sets availabl e at low redshift. For exam- 
ple, ISchawinsk i et al.1 (|2006f ) derive an empirical relation 
for a critical black-hole mass (as a function of velocity 
dispersion) above which the outflows from these black 
holes suppress star formation in their hosts. 

Although the low redshift studies benefit from large, 
high-quality surveys and detailed information, they do 
not enable us to witness the star formation truncation 
in more massive galaxies, as stellar populations studies 
show that thes e objects formed most of their stars a t high 
redshift (e.g., Ivan Dokkum fc van der M arcl 200(| and 
references therein). Recent attempts to directly witness 
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the AGN feedback process at z > 2 have been limited to 
studies in which galaxies ar e selected for their stro ng nu- 
clear activity. For example, iNesvadba et al.1 (|2006P argue 
that an AGN driven wind is the only plausible mecha- 
nism to explain the outflow of the gas seen in a powerful 
radio galaxy at z = 2.16. A more statistical approach 
is studying the X-ray properties of AGNs with redshift. 
These studies found that AGNs show a top-down behav- 
ior, such that the space de nsity of the more luminous ones 
peaks at higher redshift (ISteffen et al.|[200^ : lUeda et al.l 
l2003HHasinger et al.ll2005fl . iHeckman et all (|2004D claim 
that this behavior reflects the decline of the characteris- 
tic mass-scale of actively accreting black holes with red- 
shift. The fact that this behavior is strikingly similar 
to what is found for the stellar populations of galaxies, 
such that the stars in more massive galax i es are formed 
at hi gher redshift (e.g.. lCowie et al.lll996t Uuneau et all 
20031), could be a clue that the two are strongly related. 
In order to relate these two behaviors, and understand 
the role of AGNs in the SFH of massive galaxies it is cru- 
cial to study massive galaxy samples at earlier epochs. 

At a redshift of z ~ 2.5 massive galaxies (> 
lO^ Mp)) range from sta r bursting to evolved 12 systems 



TABLE 1 

Sample 



e.g.. iFranx et al 



Labbe et al. 2005 



2003; iForster Schreiber et all 12 004; 
_ .Reddv et alj 120061; iPapovich et al.l 
20061; IKriek et alj l2006bl ; iWuvts et"afl 120061 ). Thus a 
massive galaxy sample in this redshift range may be 
expected to contain all evolutionary stages of the pro- 
cess that transforms a massive star forming galaxy into 
a red, quiescent system. As massive galaxies are bright 
at rest-frame optical wavelengths, a representative mas- 
sive galaxy sample can be obtained by selecting at near- 
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Furthermore, for de- 
the star formation and nu- 



infrared wavelengths (e.g., 
I2004bt Ivan Dokkum et al. 
tailed information about 
clear activity in the massive galaxies, spectroscopic in- 
formation is required. As the average massive galaxy at 
2.0 < z < 3.0 is faint in t he rest-frame UV (Rab = 25.9, 
Ivan Dokkum et aL1l2006l ). it is beyond the limits of op- 
tical spectroscopy. Thus, if we want to obtain spectro- 
scopic data on a representative massive galaxy sample at 
z ~ 2.5, we need to observe at NIR wavelengths. 

In order to understand the formation of massive galax- 
ies, we are undertaking a NIR spectroscopic survey for 
massive galaxies at z ~ 2.5 with GNIRS on Gemini- 
South. Ideally we would study a mass-limited sample 
with no regard for luminosity or color. However, many 
massive galaxies are too faint for NIR spectroscopy on 
today's largest telescopes. Therefore, we study a K- 
selected sample, which is much closer to a mass-limited 
sample than an i?-selected sample. We note, however, 
that we may miss massive galaxies with comparatively 
high M/L ratios. So far, our spectroscopically confirmed 
if-selected sample consist of 20 galaxies with 2.0 < z < 
2.7. Nine of these galaxies show no emission lines and are 
characterized by strong B almer/4000 A breaks . These 
galaxies are discussed in IKriek et al.l (|2006bl ). Here 
we discuss GNIRS and follow-up VLT/SINFONI ob- 
servations of the emission line galaxies in the sample. 
Throughout the paper we assume a ACDM cosmology 
with n m = 0.3, n A = 0.7, and H Q = 70 km s" 1 Mpc" 1 . 

12 In this paper, "evolved" is shorthand for having a low specific 
star formation rate (SFR/M*) 
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a Integration times for SI NFONI. b Integratio n time s for G NIRS . c Ob- 
servations presented by Ivan Dokkum et al.1 lt2005h and Ikriek et ail 
l|2006aD . The R magnitude is adopted from lDaddi et alj d2004alh 

All broadband magnitudes are given in the Vega-based 
photometric system. 

2. DATA 
2.1. Sample 

The galaxies studied in this paper are drawn from 
a spectroscopically co nfirmed if -selected galaxy sam- 
ple at 2.0 < z < 2.7 (|Kriek et al.l l2006bl ). This spec- 
troscopic sample was originally selected from the op- 
tical and deep NIR infrared photometry provided by 
the Multi- wavelength Survey by Yale-C hile (MUSYC, 
iGawiser et al] 120061; lOuadri et al.ll2007fh The original 
sample, selected for 2.0 < z p hot < 2.7 and K < 
19.7, contains 26 galaxies all observed wit h the Gem- 
ini N ear-InfraRed Spectrograph (GNIRS, Elia s et ail 
2006). Twenty galaxies have a spectroscopic redshift 
within the targeted redshift range. Our if-selected sam- 
ple seems representative for all K < 19.7 galaxies in 
this redshift range, as both a Mann- Whitney U-test 
and a Kolmogorov-Smirnov test show that the spectro- 
scopic sample has a similar distribution of rest-frame 
U — V colors as the large mass-limited photom etric sam- 
ple (> 10 n M Q ) by Ivan Dokkum et all (|2006D when ap- 
plying the same if -magnitude cut. 

Emission lines were detected for eleven of the twenty 
galaxies. The nine galaxi es without emission lines are 
presented and discussed in Kriek et al. (2006b). We ob- 
served ten of the eleven emission line galaxies with SIN- 
FONI (jEisenhauer et all l2003t iBonnet et all l2004h . to 
obtain higher resolution spectra, and two-dimensional 
(2D) information on the line emission. For one of the 
emission line galaxies we took no follow-up SINFONI 
spectra, as this galaxy is already discussed in detail by 
Ivan Dokkum et al.l (|2005h . 

2.2. GNIRS spectra 

The original spectroscopic galaxy sample was observed 
with GNIRS in 2004 September (program GS-2004B- 
Q-38), 2005 May (program GS-2005A-Q-20), 2006 Jan- 
uary (program GS-2005B-C-12) and 2006 February (pro- 
gram GS-2006A-C-6). We used the instrument in cross- 
dispersed mode, in combination with the short wave- 
length camera, the 32 line mm" 1 grating (i? ~ 1000) 
and the 07675 by 6'.'2 slit. In this configuration we ob- 
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Stellar population properties of the emission line galaxies 
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Note. — The stellar population properties arc derived from fitting the low-resolution con- 
tinuum spectra and optical photometry by stellar population models. The errors present the 
68% confidence intervals derived using 200 Monte Carlo simulations. 



tained an instantaneous wavelength coverage of 1.0 - 2.5 
(Mm. The integration times for the emission line galaxies 
are listed in Table 1. The observational techniques and 
red uction of th e GNIR S spectra are described in detail 
by iKriek et"aT1 (|2006a) . For each galaxy we extract a 
one-dimensional (ID) original and low-resolution binned 
spectrum. 

We use the low-resolution continuum spectra in com- 
bination with the optical photometry to obtain the in- 
tegrated stellar population properties of the galaxies. 
The near-infrared spectra are flux calibrated using JHK 
ph otometry. The spe c tra a nd UBVRIz fluxes are fit 
by iBruzual fc Charlotl (|2003l ) stellar population models 
with exponentiall y declining SFHs, f ollowing the tech- 
nique de s cribed in IKriek et alj (|2006al lbl). We assumed a 
Salpeterl (119551) IMF and adopted the reddening law by 
Calzetti et al.1 (|2000f) . The redshift was set to the emis- 
sion line redshift during fitting. We allowed a grid of 24 
ages between 1 Myr and 3 Gyr, 40 values for Ay between 
and 4 mag, and 31 values for r (the SFR decaying time) 
between 10 Myr and 10 Gyr. 

GNIRS is uniquely capable of this technique as for 
z ^ 2.5 galaxies the instrument covers the whole rest- 
frame optical wavelength regime in one shot, from blue- 
wards of the Balmer break up to 7000 A. The stellar 
population modeling, driven by the optical continuum 
break, yields redshifts as well. This is in particular use- 
ful for galaxies without detected emission lines, such as 
passive systems. The stellar population properties of the 
emission line galaxies studied in this paper are listed in 
Tabled 

2.3. SINFONI Spectra 

We observed the ten emission line galaxies with near- 
infrared integral-field spectrograph SINFONI during two 
runs, on 2005 December 10-13 (076.A-0464) and 2006 
March 3-4 (076.A-0718). The weather conditions during 
both runs were fairly stable, with a median seeing of 0'.'4 
in the NIR. We use the H+K grating {R ~ 2000, A = 
14500-25000A) over the 8"x 8" field of view (FOV). The 
FOV is sliced into 32 slitlets, and the spatial sampling 
corresponding to this configuration is 0'.'25 x 0'.'125. We 
observed the galaxies according to an ABA'B' on-source 



dither pattern. The offset between A and B is half the 
FOV (4") in the directions perpendicular to the slitlets. 
The offset between A and A' is only 1". This dither 
pattern enables an accurate background subtraction, as 
we observe empty sky on both sides of the object in each 
slitlet. 

The raw SINFONI spectra were reduced using custom 
IDL scripts that perform the following steps. We start 
by correcting for the detector response by dividing by 
the lamp flats. Next, we determine the positions of the 
spectra on the detector for each slitlet from the "distor- 
tion frames" . These frames are constructed by moving 
an illuminated fiber in the direction perpendicular to the 
slitlet over the FOV. For every slitlet spectrum we de- 
termine the position with respect to the other spectra at 
each wavelength by tracing the illuminated fibre. The 
slitlet lengths, which are different for each slitlet, are de- 
rived from the obtained fiber traces in combination with 
the sky emission from the raw object frames. 

Next, we make a combined bad pixel map for each in- 
dividual frame, that identifies cosmic rays, hot pixels and 
outliers. To identify cosmic rays, we first have to remove 
the sky emission. An initial sky removal is performed 
by subtracting the average of the previous and succes- 
sive frame. The remaining sky residuals are removed by 
subtracting the median flux at each wavelength for each 
slitlet spectrum, using the derived locations of the spec- 
tra on the detector, and masking the object position. 
Next, we identify cosmic rays on the obtained images 
using L. A. Cosmic (Ivan Dokkumll200lh . We add any re- 
maining 4er outliers to the map as well. This map is 
combined with a common bad pixel map, constructed 
from flat and bias frames. In what follows the combined 
bad pixel map will be transformed in the same way as 
the science images. 

We return to the raw images and again we perform a 
simple sky-subtraction, now using the bad pixel map to 
reject cosmic-rays and other defects. Next, we cut and 
straighten all spectra, using the derived positions. For 
each slitlet spectrum, we straighten the skylines and per- 
form the wavelength calibration in one step, so that the 
data are resampled only once. Now, we can accurately 
remove remaining sky at each wavelength for each slitlet 
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Fig. 1. — One dimensional spectra in the wavelength region around H/3 and [Om] AA4959, 5007, and Ha and [Nil] AA6548, 6583 of ten K- 
selected emission line galaxies at 2.0 < z < 2.7. The wavelength is presented in rest-frame and the flux is given in 10 -19 ergs s _1 cm _2 A — 1 . 
The vertical dotted lines present the positions of the expected H/3, [Om] AA 4959, 5007, [Nil] A6548, Ha and [Nil] A6583 lines. The red 
line presents the best fit to the three emission lines. The green line is the best continuum fit to the low-resolution GNIRS spectra. Gray 
shaded areas present the noise spectrum. For several galaxies we combined the SINFONI with the GNIRS spectra for reasons explained in 
the text. 



spectrum, by masking the object spectrum. The previ- 
ous steps yield a 3D datacube for each exposure. Finally, 
we combine the data cubes of the individual exposures, 
using the bad pixel cubes and the offsets. The final cube 
is divided by a response spectrum, created from the spec- 
tra of AVO stars, matched in airmass, and reduced in a 
similar way as the science objects. 

2.4. Extraction of one- dimensional spectra 

We extract the ID spectra in two different ways de- 
pending on whether the Ha and [N n] emission lines are 
blended due to a strong velocity gradient. For these 
galaxies (ECDFS-3694 and ECDFS-10525) we remove 
the relative velocity shifts along the spatial direction be- 
fore constructing a ID spectrum, in order to obtain more 
accurate line measurements and ratios. While informa- 
tion on the total velocity dispersion is thereby lost, this 
procedure is justified for the purpose of extracting the 
integrated line fluxes for the analysis presented in this 
paper. 

For both methods we start by binning the final SIN- 
FONI cubes in the spatial direction, to combine the 
'spaxels' of 0"I25 by 0"25 to spatial elements of 0"25 by 



0'.'25. As we are mainly interested in studying the spec- 
trum of the line-emitting gas rather than the continuum 
emission, we optimize the extraction of the ID spectra 
in the wavelength range around Ha and [N n] . We make 
an average reconstructed image of the data cube in the 
wavelength region around these lines. The wavelength 
region within twice the velocity dispersion of the lines 
is included, avoiding wavelengths with strong OH lines 
or low atmospheric transmission. Note that the velocity 
dispersion is still unknown at this stage, so we have to it- 
erate a few times to obtain the correct extraction region 
and the velocity dispersion. For galaxies with a velocity 
gradient, we take into account the relative velocity shift 
as well to determine which wavelength region is included 
Next, we select all adjacent pixels in the reconstructed 
image with a flux exceeding 0.20-0.55 times the flux of 
the spatial element with the most signal. This threshold 
value is dependent on the S/N of the line emission in 
the reconstructed images. In case no velocity gradient is 
present, the ID spectra of the selected spatial elements 
are combined to form the final ID spectrum. A ID noise 
spectrum is constructed from the spectra of all spatial 
elements that do not exceed this threshold. 
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Fig. 2. — 3" by 3" images (~25 by 25 kpc) of the continuum (color scale) and line (contours) emission. The continuum emission is 
derived from the median collapsed H+K SINFONI cubes, excluding wavelength regions with low atmospheric transmission or strong sky 
emission. The linemaps include both the Ho and [N II] A6583 emission. The contours represent the 50% and 90% of the maximum line 
emission in the galaxy. Three of the four AGN candidates exhibit compact line and continuum emission. 

For the galaxies with a velocity gradient, we fit the 
spectra of all selected spatial elements following the fit- 
ting procedure described in § 12.51 To avoid the residuals 
of skylines as being interpreted as lines, and to reduce the 
number of degrees of freedom, we fix the width to the 
velocity dispersion of the final spectrum. The allowed 
redshift is also constrained to a certain range, depend- 
ing on the maximum velocity shift of the galaxy. Note 
that these latter two requirements need a few iterations. 
We measure Ha and the two [N n] lines in the spectra of 
all selected spatial elements, and retain the elements for 
which Ha and [Nil] A 6583 have a total signal-to-noise 
(S/N) > 5. We visually inspect the spectra of each ele- 
ment to check if the lines are correctly interpreted. The 
rest-frame ID spectra of these selected spatial elements 
arc averaged together to form the final ID spectrum. 
We constructed a final noise spectrum by quadratically 
adding the rest-frame noise spectra. 

Finally, we added the ID spectra to the previously 
measured GNIRS spectra in cases where the S/N is low 
and the SINFONI data show no evidence for strong ve- 
locity gradients. 

2.5. Line measurements 

We obtain emission line ratios, velocity dispersions and 
equivalent widths of the emission lines by modelling the 
extracted ID spectra. The H and K spectral bands were 
fit separately. The H band contains the lines H/3 and 
[Om] AA4959, 5007. The K band covers Ha and [Nil] 
AA6548, 6583. We fit Gaussian models to each set of lines 
simultaneously, assuming a similar width for all three 
lines and one redshift. We note that the different emis- 
sion lines may not have similar line width as they can 
originate from different physical processes. However, as 
we will use the emission line measurements to identify 
the main contributor to the line emission, we assume 
that all lines originate from the same emission line re- 
gion. Furthermore, the [Nil] and Ha lines have similar 
widths for the galaxies for which we could measure the 
lines separately. 

We adopt the ratios of transition probablity between 



the two [Nil] lines and the two [Om] lines of 0.34 and 
0.33 respectively. Thus, for each fit there are four free 
parameters: redshift, line width, and the fluxes of the 
two emission lines. The line widths were fixed to the 
best-fit value obtained from the set of emission lines (ei- 
ther H/3 and [Om] AA4959, 500 or Ha and [Nn] AA6548, 
6583) with the highest S/N. For the continuum we use 
the best fit to the GNIRS low-resolution continuum spec- 
tra, corrected for velocity broadening, assuming that the 
stellar dispersion is similar to the gas dispersion. Thus 
the continuum directly includes the Balmer absorption. 

The ID spectra are fitted by minimizing the absolute 
residuals from the fit weighted by the noise spectrum. 
This fitting method is preferred over \ 2 fitting, as it min- 
imizes the influence of sky lines and strong noise peaks. 
Errors on the flux measurements were determined by fit- 
ting 500 simulated spectra, that were constructed from 
the original spectrum and the photon noise. In these 
simulations we also varied the continuum according to 
the probability distribution that followed from modeling 
the GNIRS spectra. In cases where we fix the width, as 
derived from a brighter emission line in the same spec- 
trum, we also vary the assumed width in the simulations 
according to the corresponding probability distribution. 

For all galaxies except 1030-2026, we start by fitting 
Ha and the [Nil] lines. We use the obtained redshift, line 
width and its probability distribution to fit H/3 and the 
[O III] lines. For 1030-2026 Ha is strongly blended with 
the [Nil] lines, even in the higher resolution SINFONI 
spectra. As this galaxy has a clear [O ill] A5007 detection 
in both SINFONI and GNIRS spectra, we swap the or- 
der. Thus for this galaxy we first derive the redshift, the 
line width, and its probability distribution from model- 
ing [Om] and H/3. Using the derived redshift, line width 
and its probability distribution, we can now measure Ha 
and the [Nil] lines. 

Figure [1] shows the ID spectra and best-fit model for 
all detected emission lines. The final ID spectra include 
the SINFONI data, or a combination of both SINFONI 
and GNIRS. As the SINFONI extraction method is op- 
timized for the line emitting gas, and the GNIRS for 
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TABLE 3 
Emission line modeling results 



id 


<7 a 


W Hoi h 


W[NII] 




Wfom] 


[Nn] /Ha 


[O m]/H/3 c 


[Om]/H/3 d 


1030-807 


115+21 

—52 


61.7+JV 

— 8.6 


18.7+?^ s 

— 14.5 


13.0+ 6 ,* 

— 5.9 


— 2.0 


0.33+ n '°,5 
— 0.25 


< 0.83 


< 0.86 


1030-1531 


79+54 


62.8l 4 8 9 8 


7 6+ 61 
'■"-12.2 


9.51JJ 




< 0.35 


1 o 7 +0.31 

-0.37 


< 1.99 


luOU-ZUZD 


404+117 
434 _ 64 


d2 " y -12.5 


91 n+ n 1 


9 n~^3.2 

2 - U -5.2 


14 .+ 4 ' 6 
— 1.1 


n R4+ - 61 

U - b4 -0.20 


^> 1 .O^i 


S Q 70 


1030-2329 


80tf 3 


51.1±«;2 




7 fi+ 2 ' 8 
— 1.4 


9+5.2 
°- -2.5 


71+ 08 
u -' 1 -0.06 


< 2.30 


< 2.51 


1030-2728 


114+29 


32.4+»-» 


19.4+J-* 




Q+ 2 ' 5 


63+° 17 




< 1.77 


ECDFS-3662 


141±|g 


54.6±iV 


29 - 4+ 68 




1 s+°- 4 


5fi+ 09 




< 1.06 


ECDFS-3694 


139±g 


95.8+f 7 5 


41.511° 


28.7+iV 


48.9l 3 ° 4 7 


45+O.05 


9 (11 +1-28 


< 4.88 


ECDFS-3896 




31.5±&* 




r 7 +1.9 
°-'-5.0 


9 9+1.9 
z - -2.3 


1 09+ ' 31 




< 1.32 


ECDFS-5754 


160±JJ 


76.9l 7 -^ 




99 q + 9-8 
zz '°-11.7 


°- i -6.0 


20+° 19 


< 1.87 


< 2.21 


ECDFS-10525 


202+29 




57.9l?'I 


7.8l^ 2 o 7 


-5.21^ 


™+°- 06 




< 0.35 


CDFS-6036 C 




99+10 

»»_io 


601^ 


12±! 


99lJo 


n co+o.09 


8-28111 





Note. — The errors present the best 68% confidence intervals derived using 500 Monte Carlo simulations. All upper 
and lower limits are la. The emission lines ratios present the flux ratios. 

a In units of kins' -1 , corrected for instrumental resolution. For ECDFS-3694 and ECDFS-10525 the velocity gradient is 
removed from the velocity dispersion, to deblend the emission lines. All equivalent widths are in rest-frame, corrected for 
Balmcr absorption, and given in A. c Derived from the emission line modeling.^ Derived from the lower limit on H/3 and 
the upper limit on [Olll]. The lower limit on H/3 is estimated using the lower limit on Ha, the intrinsic ratio between Ha 
and H/3 and the modeled upper limit on the continuum attenuation. Furthermore we assume extra attenuation towards 
H II regions. 6 Measured by Ivan Dokkum et alj J2005D 



the continuum, the final ID spectra may not be repre- 
sentative for the whole galaxy. Nevertheless, a compari- 
son between GNIRS and SINFONI for the galaxies which 
could be measured separately, show that the Wn a 's are in 
good agreement: ECDFS-10525: Whq.gnirs = 1161^ A 
and W Hq ,sinfoni = 117l 4 3 A; 1030-807: W Ha ,GNiRS = 
56±? 3 A and Wh^sinfoni = 62^ 3 A; ECDFS-3694: 

WWgnirs = 125l^ 5 A and W Hq , S infoni = 96^ 7 A. 

The best-fit equivalent widths, line widths and emis- 
sion line ratios are listed in Table [3] The values for the 
emission line ratios are given when both lines have a 

> 2(T detection. In case only one of the two lines has a 

> 2a detection we give a 2a upper or lower limit. As 
can be seen in Figure [1] Ha is detected for all galaxies, 
and [N n] A6583 can be measured for nine of ten galax- 
ies. Thus, we can determine the value for [Nil] /Ha for 
these nine galaxies. For 1030-1531, which has no [Nil] 
detection we give a 2a upper limit. For only two galax- 
ies both H/3 and [O in] A5007 are detected at > 2a. For 
these galaxies the ratio [O III] /H/3 can be measured di- 
rectly. For four objects we obtained lower or upper limits 
from the fitting procedure, as one of the two lines (H/3 
or A5007) had a 2a detection. For example, for galaxy 
1030-807 the H/3 line is detected at 3 sigma and [O in] 
A5007 has an upper limit. Thus, we can derive a 2a up- 
per limit on [O lil]/H/3, from the 97.5% cut of the best-fit 
[O m] /H/3 values of the 500 simulations. For the remain- 
ing four galaxies the modeling results yielded no limit on 
[O in] /H/3, as neither of the two lines was detected. 

For all galaxies we apply a second method to constrain 
[Olll] A5007/H/3, based on the intrinsic ratio of Ha/H/3 
(= 2.76). We determine a 2a lower limit on H/3 from the 
2a lower limit on Ha. Furthermore, we attenuate H/3 
using the 2a upper limit on the best-fit modeled Ay and 
assuming extr a extinct i on towards Hii regions ( factor 
of 0.44 in A v , lCalzettilH997l : iCalzetti et all [20001. We 
combine the attenuated 2a lower limit on H/3 with the 
modeled 2a upper limit on [O in] A5007 to derive the 2a 
upper limit on [O in] /H/3. The limits are listed for each 



galaxy in the last column of Table [3l 

The Ha line of 1030-2026 is not well fitted in Figure [TJ 
This is not caused by a wrong redshift measurement, as 
we know the redshift of this galaxy very accurately from 
the [O in] A5007 line. There seems to be a second peak at 
a rest-frame wavelength of 6600 A. Although this "line" 
falls on top of a strong OH line, the emission feature looks 
real in the 3D cube. It could well be Ha emission from 
a companion galaxy or a star forming region located in 
the outer parts of the galaxy. 

2.6. Line maps 

For each galaxy we make reconstructed images of the 
continuum and line emission separately. The continuum 
images includes all observed H+K emission, excluding 
wavelengths with low atmospheric transmission or strong 
sky emission. The linemaps include both the Ha and the 
[Nil] A6583 emission. The S/N is not high enough to 
make linemaps of these two components separately. We 
included the line emission within two times the velocity 
dispersion of the emission lines, accounting for a velocity 
gradient if present. The continuum emission is removed 
using the flux measurements in the wavelength range sur- 
rounding the lines. Finally, we convolve the line map by 
a boxcar of 3 pixels (0'.'375), to smooth out the noise. 

Maps of the continuum and line emission are presented 
in Figure [2] The spatial sampling of the reconstructed 
images is 0'.'125 by 0'.'25, unlike the reconstructed images 
used to extract the spectra, which are binned to 0'.'25 by 
0'.'25. The area over which the ID spectra are extracted is 
about similar to what is included by the outer contours. 

3. ORIGIN OF THE LINE EMISSION 

A key goal of this paper is to elucidate the origin of the 
Ha emission detected in slightly over half of the massive 
galaxies in our sample. In the local universe, Ha emis- 
sion is usually an indicator of star formation, and the 
luminosity of Ha directly correlates with t he instanta- 
neous star formation rate (|Kennicuttl Il998f ). However, 
other processes can also ionize Hydrogen, in particular 
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Fig. 3. — Diagno stic diagram for spectral classification of AGNs and star-forming galaxies. The gray-scale presents the locus of ~400 
000 SDSS galaxies ( Kauffmann et al.l 2003b; Trcmonti ct al. 200J). The orange solid lin e is the empirical division between galaxies for 
which the line emission originates from Hn regions and AGNs for the SDSS galaxies by Kauffmann ct al. (2003b). The orange dashed 
line presents the theoretical upper limit by Kewley et all (12001) for star forming galaxies. Galaxies between these two dividing curves are 
classified as composite H II- AGN gal axies bv lKewlev et al.l (120061 ). The orange dotted line presents the division between Seyfert 2s and 
LINERs by Kauffmann ct al. (2003b). The red filled squares present galaxies with a specific SFR (derived from modeling the continuum 
spectra) less than 0.05 Gyr" 1 , and the blue filled circles galaxies with higher specific SFRs (> 0.05 Gyr~ x ). Furthermore, the green squares 
indicate galaxies with compact line emission. The galaxies with black circles are identified as AGN candidates, based on their [Nil] /Ha 
ratios, spatial extent of the line emission, and ancillary data. Further details are in the text. All upper limits are 2<r and the error-bars are 
all lcr. For 1030-2026 we have both a 2cr upper and lower limit. 



hard radiation from a central active nucleus. In this sec- 
tion we use different indicators for star formation and 
AGN activity to determine the nature of the population 
of emission-line galaxies. 

3.1. Emission line ratios 

Local star-forming galaxies in the SDSS follow well- 
defined tracks in diagnostic diag rams featuring var- 
ious emission line ratios (e.g. iBaldwin et al.l [1981L 
iVeilleux k. Osterbro ck 1987). For our sample, the appro- 
priate diagram is [Oin] A5007/H/3 versus [Nil] A6583/Ha 
as these lines are relatively strong and are covered in 
both our GNIRS and our SINFONI spectra. This dia- 



gram is shown in Figure 3. Galaxies for which lines origi- 
nate from photo-ionization by young stars in H II regions 
fall on the well-defined sequence. Local galaxies outside 
this metallicity driven sequence are dominated by other 
ionization sources, typically photo-ionization by a hard 
s pectrum such as p r oduced by an AGN. 

Kauffma nn et al.l (|200 3h) empirically separates the 
AGNs from the H n sequence by the solid line in Figured 
The extreme starb urst classification line derived by 
iKewlev et al.l ([2001) is presented by the dashed line. The 
galaxies between these two dividi ng lines are classifie d 
as composite Hi i-AGN galaxies bvlKewlev et al.l f2 006'l. 
However, both lErb et all (|2006af ) and iShaplev et al.l 
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Fig. 4. — The measured Wn a vs. the expected Wn a derived from 
the sp ectral continuum shape (see Figure [5} and the IKennicut3 
(1998) relation between SFR and Ho luminosity. In case Hn re- 
gions are the only contributors to the line emission, a galaxy falls 
on the expected 1-to-l relation (dotted line). Additional extinction 
towards star forming regions moves a galaxy downwards of the re- 
lation. Objects that fall above the dotted line most likely have 
another ionization source contributing to the line emission as well. 
The symbols are similar as in Figure [3] 



(2005) suggest that the Hn driven sequence is offset to 
the right at high redshift. Thus, we have to be careful 
when applying the classification scheme as derived for lo- 
cal galaxies to our high-redshift sample, as the behavior 
of the ionization ratios at high redshift is not well under- 
stood. The galaxies outside the Hn sequence are gen- 
erally divided into LINERs and Seyfe rt 2s. The dotted 
line in Figure [3] is the dividing line by [Kauffmann et all 
(|2003bf) . 

Our eleven emission line galaxies are indicated by the 
red filled squares (specific SFR < 0.05 Gyr -1 ) and blue 
filled circles (specific SFR > 0.05 Gyr -1 ) in Figure OS 
Remarkably, only four objects fall in the re gion of galax- 
ies wi th pure star formation as defined bv lKewlev et al.l 

(2006) , while the other seven galaxies are classified as 
AGN or composite Hn-AGN. Three of these galaxies fall 
in the Seyfert 2 regime, three galaxies may be LINERs 
or Seyfert 2s, and one is classified as a LINER. 

While Seyfert 2s are generally accepted as AGNs, the 
power source of LINERs is still debated. Although 
photo-ionization by an AGN is often the most straight- 
forward explanation, LINER emission has also been ob- 
served in extra-nuclear regions a ssociated with large- 
scale outflows and relat ed shocks ()Dopita fc Sutherland! 
Il995t iLipari et al.ll2004|). Shock-producing winds can be 
drive n by AGNs (e.g.. lCecil et al.ll2000t iNesvadba et~aT 



2004). 



2006), but also by strong starbursts (e.g.. ILipari et al 



3.2. Other diagnostics of AGN activity 

Given the ambiguity in the interpretation of LINER 
spectra it is necessary to examine other indicators of 
star formation and AGN activity, to determine for which 
of the seven galaxies with non Hn region- like ratios 



an AGN is the most likely explanation (either directly 
through photo-ionization or through shocks produced by 
an AGN-driven outflow). We first consider which are the 
most important additional diagnostics at our disposal: 

X-ray emission: As is well known AGNs can be ef- 
ficiently identified by their X-ray emission, which is 
thought to be due to up-scattered UV photons from 
the accretion disk. AGN-induced X-ray emission can 
be distinguished from that induced by star formation 
by the hardness ratio and (particularly) the luminos- 
ity. For the galaxies in SDSS1030 we use XMM data 
with a depth of 100 ks (Uchiyama et al. 2006, in prepa- 
ration). Three of the ECDFS galaxies (ECDFS-3662, 
ECDFS-3694 and CD FS-6036) are in the CDFS proper 
(jGiacconi et al.ll2002f ) , for which very deep 1 Ms Chandra 
data are available. The other CDFS galaxies (ECDFS- 
3896, ECDFS-5754 and ECDFS-10525) are in the "Ex- 
tended" CDFS, fo r which we use 250 ks Chandra data 
(|Virani et alJ 12006). Interestingly, only one out of seven 
galaxies is detected. Limits are given in Table 1. These 
limits indicate that the AGNs, if present, are either 
highly obscured or accrete at sub-Eddington rates. 

Compactness: The spatial distribution of the line emis- 
sion can provide information on its origin, as narrow line 
regions of AGNs are generally compact. However, ex- 
tended line emission does not rule out the presence of 
an AGN for galaxies with active star formation, as star 
forming regions contribute to the line emission. Further- 
more, AGNs can produce outflows, which may result in 
extended line emission. As can be seen in Fig. [2j in sev- 
eral galaxies the line emission is very compact whereas 
in others it is extended. 

Ha equivalent width: As we have estimates for the SFR 
in these galaxies from the stellar continuum fitting, we 
can compare the observed strength of Ha to that pre- 
dicted from the best-fit synthetic spectrum. An excess 
of Ha emission could indicate another ionization source 
than star formation, such as AGN activity. We com- 
pute the predicted Ha line luminosity using the SFR 
from the best-fit model to the observed continuum emis- 
sion (which inclu des Ay as a free parameter) and the 
iKennicuttl (|1998| ) relation. We then divide this estimate 
by the continuum luminosity density around the wave- 
length of Ha as determined from the best-fit synthetic 
spectrum and corrected for the best-fit extinction 

In Figure [4] we compare the measured and expected 
Who- The relation between the two is dependent on the 
dust geometry, such that extra extinction towards H II 
regions would move a galaxy below the 1-to-l relation 
(which is indicated by the dotted line). As it is unlikely 
that the attanuation in sta r forming sites is less than t he 
continuum attenuation fsee lCid Fernandes et a l. 2005), a 
galaxy is not expected to lie above the dotted 1-to-l rela- 
tion. However, in case of another ionization source such 
as AGN activity, the measured Wh q would be higher 
than the expected. We indeed find that several galax- 
ies fall above the expected relation, which implies that 
an ionization source other than star formation may be 
contributing to line emission as well. 

Star formation activity: If a galaxy has a high SFR the 
current starburst might produce an outflow which results 
in AGN-like emission line ratios. Although this outflow 
is hard to confirm, the absence of a starburst immedi- 
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ately rules out this scenario, and infers that ionization 
ratios are indeed due to AGN activity. To indicate which 
galaxies have an engine for such an outflow, we divided 
the sample in galaxies with high (blue filled circles) and 
low SFRs (red filed squares) in Figure 3. 

Next, we assess for each of the seven candidates 
whether the line emission is most likely dominated by 
an AGN or some other process, using these criteria and 
others particular to individual objects. We stress that 
in none of the cases it is completely clear-cut: even in 
the local universe, with the availability of vastly superior 
data to ours, it is often impossible to cleanly identify the 
relative contributio ns of AGNs and star formation to th e 
line emission (e.g.. lFilho et al.ll2004 fKewlev et alJl2006t ). 

1030-2026: This galaxy has Seyfert 2 emission line ra- 
tios, compact line emission, and a low SFR (as implied by 
the continuum modeling). Thus, a starburst driven wind 
is very unlikely for this galaxy. Furthermore, Figure [4] 
shows that ongoing star formation, as derived from the 
stellar continuum, cannot account for the observed Ha 
emission. Finally, the velocity dispersion of ~450 kms -1 
of [O in] A5007 may be indicative of an AGN. We classify 
this galaxy as an AGN. 

1030-2329: The emission line ratios of this galaxy are 
indicative of an AGN, or a composite H n-AGN galaxy. 
The excess of Ha emission in Figure [4] supports the pres- 
ence of another ionization source. Furthermore, for this 
galaxy a starburst driven wind is very unlikely, due to the 
combination of a low un-obscured SFR (as implied by the 
continuum) and compact line emission. We include this 
galaxy in our AGN selection. 

1030-2728: This galaxy has similar supporting diag- 
nostics as 1030-2329, but the line emission is not com- 
pact. However, the spectrum shows hints that the [Nil] 
emission peaks in the central emission blob in Figure [21 
while the lower-right emission peak is dominated by Ha 
emission. Unfortunately, the S/N of the individual line 
maps is not high enough to draw firm conclusions. As 
all other diagnostics indicate an AGN as the dominat- 
ing ionization source, we add this galaxy to our AGN 
selection. 

ECDFS-3662: This galaxy probably falls in the com- 
posite H II- AGN region of Figure El Unfortunately, we 
cannot discriminate between shock-ionization by a star- 
burst driven outflow and an AGN, as this galaxy has a 
high SFR and extended line emission. Other diagnostics 
do not provide evidence for an AGN: the object has no X- 
ray counterpart (in the 1 M s data), and is not identified 
as an AGN (nor ULIRG) bv lAlonso-Herrero et all (|2006h 
using the infrared continuum shape. Furthermore, the 
observed Ha emission is consistent with star formation. 
Although we cannot rule out the presence of an AGN, 
we will not include this galaxy in our AGN sample. 

ECDFS-3694'- This galaxy is classified as a composite 
H II- AGN galaxy in Figure [31 It has a high SFR, ex- 
tended line emission, and shows a large velocity gradient 
of ~1400 kms -1 . This galaxy falls in the field examined 
bv lAlonso-Herrero et al.l (|2006l ). and is not classified as a 
ULIRG or AGN according to its mid-infrared SED shape. 
It is undetected in the 1 Ms X-ray imaging. A starburst 
driven outflow seems the most plausible scenario for this 



galaxy, although an AGN cannot be ruled out. Thus, this 
galaxy will not be part of our AGN candidate sample. 

ECDFS-3896: The emission line ratios for this galaxy 
are indicative of an AGN or a composite Hn-AGN 
galaxy. This galaxy has a high SFR derived from the 
stellar continuum, which is consistent with the observed 
Ha emission. Nevertheless, a starburst driven wind is 
somewhat unlikely, as the line emission shows a very com- 
pact structure and the galaxy shows no velocity gradient. 
Furthermore, as will be discussed in the next section, the 
UV-emission which boosts the SFR in the model fits, may 
well be due to continuum emission by the AGN. Optical 
spectroscopy is needed to clarify this situation. As an 
AGN seems the most plausible cause for the emission 
line ratios, we include this galaxy in our AGN sample. 

CDFS-6036: This galaxy seems t he most convincing 
Seyfer t 2 in Figure [3j Nevertheless, Ivan Dokkum et all 
(2005), who studied this object in detail, suggest that the 
ratios are caused by shock ionization due to a starburst- 
driven wind. The main evidence is the extension of the 
high [N n] /Ha ratios to the outer parts of the galaxy and 
the presence of a strong starburst to drive the outflow. 
Furthermore, the galaxy shows a strong velocity field, 
and both the Ha line emission and the X-ray detection 
in the soft band (see Table 1) are consistent with the 
SFR as derived from the stellar continuum (see also Fig- 
ure [3}. The rest-frame UV spectru m of this galaxy doe s 
not show any AGN features either ()Daddi et alj |2004a). 
The galaxy has a power -law SED (a = —1.2) , indic ative 
of an ULIRG or AGN (jAlonso-Herrero et al J 120061) . As 
the available diagnostics are compatible with a starburst- 
driven wind, we conservatively do not include this galaxy 
in the AGN candidate sample. We note that including 
this object in the AGN sample would not change our 
conclusions. 

3.3. Summary and AGN Fraction 

We find that of the eleven emission line galaxies only 
four have emission line ratios consistent with th e SDSS 
star forming sequence ([Kauffmann et all l2003bT ) . The 
immediate implication is that other processes than "nor- 
mal" star formation in H II regions play a major role in 
massive galaxies at z ~ 2.5. Determining the nature 
of these processes is difficult, as many different physical 
mechanisms can produce very similar line diagnostics. 
Using a variety of indicators we argue that for four of 
the galaxies the most likely source of the Ha emission is 
an AGN. These galaxies are indicated by large solid cir- 
cles in Figures 3 and 4, and are labeled "AGN" in Figure 
2 as well. They are all narrow line AGNs, as the veloc- 
ity dispersions of the emission lines are less than 2000 
kms" 1 (Table 3). 

According to this classification, the AGN fraction 
among our total sample of if-selected galaxies at 2.0 < 
z < 2.7 is 20% (4/20). Due to the various caveats 
in the AGN classification this fraction may be under- 
estimated. First, AGNs are easier to identify in qui- 
escent systems, than in actively star-forming galaxies, 
due to the strong contamination of line emission by H II 
regions. This not only complicates the emission line 
diagnostics, but also weakens the argument of the ex- 
cess of Ha emission compared to the stellar continua, 
as this excess is easier to detect in galaxies with low 



10 



Kriek et al. 



SFRs. Optical spectroscopy may help to uncover AGNs 
in star-forming galaxies, by identifying rest-frame UV 
high-ionization emission lines such as CIV and Hell. 
But a lso the MIR continuum shape and spectra l features 
(e.g., IStern et al.1 120061 : IWeedman et all l2006h or very 
deep SINFONI spectra to revea l the spatial distribution 
of [N III /Ha (iGenzel et al.ll2006t IF orster Schreiber et all 
2006b; iNesvadba et all 120061) can be used to identify 
AGNs in these galaxies. Furthermore, AGNs may have 
been missed if they are strongly obscured (mid-infrared 
data is available for only a few galaxies), or are just too 
faint to be detected. On the other hand, we cannot rule 
out shocks produced by a starburst driven wind as the 
origin of the line emission in at least one of the four can- 
didates; 

iPapovich etaH (|2006f ) find an AGN fraction of 
25% among a sample of distant red galaxies (DRGs, 
iFranx et alj |2"003) at 1.5 < z < 3.0, using X-ray imaging 
and the shape of the mid-infrared continuum. This result 
is consistent with our fraction, given the fact that DRGs 
make up t o 70% of the massive ga laxy population at 
2 < z < 3 ()van Dokkum et al"1l2006l). Our r esult is also 
consistent with the result of Ivan Dokkuml p00l . who 
find at least two AGN s among a sample of six emission 
line DRGs. Although iRubin et alj {2005) find an AGN 
fraction of onl y 5% among the 40 DRGs in the FIRES 
MS 1054 field (jForster Schreiber et al.ll2006af ). their cri- 
terion of Lx > 1.2 x 10 43 ergss -1 would almost certainly 
not pick up any of the four AGNs in our study. 

Our fraction of AGNs among if-selected galaxies is 
significant l y high er than among UV-selected galaxies. 
lErb et all (|2006bt ) identify AGNs in 5 out of 114 UV- 
selected galaxies, based on the presence of broad and/or 
high ionization emission lines in the rest-frame UV spec- 
tra, broad Ha li nes, or very high [Nil] /Ha ratios (see 
lErb et all |2006bD. Th i s frac tion is similar to the 5% 
found by IReddv et all (|2005h among the full spectro- 
scopic sample of UV-selected galaxies using direct de- 
tections in the 2 -Ms Chandra Deep F ield North images, 
and 3% found by Stei del e t al. (2002) among a sample of 
1000 LBGs at z ~ 3, using rest-frame UV spectroscopy. 
The comparison is complicated as all studies use differ- 
ent selection criteria. However, this cannot explain the 
substantial difference in AGN fraction between the UV 
and iVselected galaxy samples. We will return to this 
issue in § 4.2. 

4. IMPLICATIONS 

In the previous section we identified four AGN can- 
didates among the eleven emission line galaxies in our 
if-selected sample. In this section we discuss the nature 
of the host galaxies of these AGNs, and what this im- 
plies for our understanding of the role of AGNs in the 
star formation history of galaxies. 

4.1. Stellar populations of AGN host galaxies 

The low resolution GNIRS spectra and broadband 
SEDs for all four AGN candidates are presented in Fig- 
ure El The strong optical breaks for three out of four 
galaxies imply that the continuum emission in these 
galaxies is dominated by stellar light. This is less clear 
for ECDFS-3896, as the spectrum of this galaxy shows 
only a weak break. The best-fit stellar population mod- 
els to the spectrum and the optical photometry, when 




0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 



Fig. 5. — Binned "low resolution" spectra (filled squares) and 
optical- to-NIR photometry (open circles) of the 4 AGN host galax- 
ies. The best-fit stellar population models are drawn in gray. The 
detection of strong Balmer and/or 4000 A breaks in three of the 
four galaxies implies that the continuum emission in these galaxies 
is dominated by stellar light. The SED of ECDFS-3896 indicates 
active star formation in this galaxy, while the other three galaxies 
are best fitted by evolved stellar populations. 

assuming that the continuum emission originates from 
stars only, are also shown in Figure [5] The corresponding 
population properties are listed in Table [2j The spectral 
continuum of three galaxies (1030-2026, 1030-2329 and 
1030-2728) are best fit by evolved stellar population mod- 
els with low SFRs, while ECDFS-3896 is actively form- 
ing stars. However, as the continuum emission from the 
AGN might contribute significantly, the derived stellar 
mass and population properties for this galaxy are quite 
uncertain. 

The median absolute and specific SFRs (SFR per unit 
mass) of the AGN host galaxies are 9 Mq yr _1 and 0.04 
Gyr -1 respectively. To examine how the stellar popula- 
tions of the AGN hosts compare to those in other galax- 
ies in this redshift range, we divide the total if -selected 
sample into three classes : the quiescent gala xies without 
detected emission lines (jKriek et al.ll2006bD . the AGNs, 
and the remaining emission line galaxies. In Figure [6] 
we show the stacked spectra and composite broadband 
SEDs of these three classes. The stacked spectra are 
constructed from the individual GNIRS spectra, normal- 
ized over the continua between a rest-frame wavelength 
of 4000 and 6000 A. We applied a noise weighted stack- 
ing method to avoid the influence of sky lines and wave- 
length regions with bad atmospheric transmission. The 
broadband SEDs are normalized at a rest-frame wave- 
length of 5000 A. We also show the average of the indi- 
vidual best-fit stellar population models to the spectra 
and broadband SEDs in the lower panels. Note that 
for the quiescent galaxies the redshifts are derived from 
the continuum shape. As we have less precise redshift 
measurements for these galaxies, spectral features will 
be smoothed out. 

Figure [5] illustrates that the average spectrum and 
broadband SED of the AGN hosts is intermediate be- 
tween the star- forming and quiescent galaxies, although 
more similar to the latter. However, this result may be 
biased, as AGNs are easier to identify in galaxies with 



AGN in Massive z ~ 2.3 Galaxies 



11 



4 



x 3 : 




2000 4000 6000 2000 4000 6000 2000 4000 6000 

(A) (A) 



Fig. 6. — Stacked GNIRS spectra and composite broadband SEDs of the three different classes of galaxies in our sample: the nine quiescent 
galaxies without detected emission lines, the four AGN candidates, and the remaining seven emission line galaxies. The individual GNIRS 
spectra were normalized in the wavelength interval between 4000 A and 6000 A, and the individual broadband SEDs were set to unity 
at a wavelength of 5000 A. The dotted vertical lines present from left to right the location of [Oil] A3727, Ca(H), Ca(K), H5, H7, H/3, 
[Oio] A4959, [Oiii] A5007, Ha [Nil] A6583. The gray lines in the lower panels show the average of the individual best-fit stellar population 
models to the spectra and broadband photometry for each class. 



low SFRs (see ^ 13 . 3[) . Thus, we cannot draw firm conclu- 
sions on the stellar populations of the full class of AGN 
hosts. Nevertheless, the low SFRs in three of the AGN 
host galaxies may suggest that the AGN activity is re- 
lated to the suppression of the star formation. 

4.2. Stellar masses of AGN host galaxies 

In order to asses why the AGN fract ion among the UV- 
selected galaxies bv lErb et alJ (|2006a) is much lower than 
among our A"-selected sample, we plot the R — K and the 
specific SFR versus the stellar mass for both samples in 
Figure [7J indicating the AGNs in black. Remarkably, 
the stellar masses of the AGN host galaxies in both sam- 
ples are similar: 2.8xlO n M and 3.OxlO n M for the 
UV- and Af-selected galaxies respectively (for a Salpeter 
IMF). As the galaxies in the AT-selected sample all have 
similar stellar masses, an under- or overestimation of the 
AGN fraction will not change this result significantly. 
This preference for massive galaxies may be the main 
reason why AGNs are more common among A-selected 
than among UV-selected galaxies. 

The apparent mass dependence of AGN host galaxies is 
affected by selection effects. The difference in R—K color 
and specific SFRs among the A-selected and UV-selected 
AGN host galaxies is most likely caused by the different 
sample selections, as three of our AGN hosts are too faint 
in the res t-frame UV to be picked up by the BM/BX 
technique (|Steidel et al.ll2004D . Three out of seven AGNs 
in Figure have R — K > 5.2. Due to the limitations 
of both surveys, we have no information about galaxies 
with R— K > 5.2 and stellar masses below 10 Mq, and 
it could be that AGNs reside in these galaxies as well. 
However, even if we ignore all galaxies with R — K > 5.2 
a strong mass-dependency is apparent in Figure [7^,. 

One can argue that this stellar mass dependency of 



AGN host galaxies simply reflects the decrease of black 
hole mass (Mbh) when going to lower mass galaxies, 
which would make the AGNs more difficult to detect rel- 
ative to the contribution from star formation. We cannot 
address this bias fully with currently available data, but 
we can roughly estimate the effect by normalizing the 
star formation rate by the stellar mass (Fig. [7)d) . Con- 
sidering a line of constant specific SFR, and assuming 
that the black hole mass scales with the stellar mass of 
a galaxy (e.g.. iDesroches et al.l 12006). the line emission 
originating from accretion onto the black hole should be 
equally easy to detect as the line emission from star for- 
mation. However, at a specific SFR between 0.1 Gyr -1 
and 1.6 Gyr -1 we only find AGNs in very massive galax- 
ies. This implies that the black hole accretion at high 
redshift is most effective at higher masses. 

This ignores the effect of variations in extinction, which 
are known to occur in these samples. In case star forma- 
tion is more extincted at the high mass end, AGN emis- 
sion might be easier to detect in these galaxies. We also 
note that the samples are incomplete in the shaded areas 
indicated in Figure [Jj Furthermore, the stellar masses 
and specific SFR derived for ECDFS-3896 and the UV- 
selected galaxies might be affected by continuum emis- 
sion from the AGN. 

4.3. Downsizing of AGN host galaxies 

Deep X-ray studies suggest that the AGN population 
exhibits cosmic downsizing, as the space density of AGNs 
with low X-ray luminosities peaks at lower redshift than 
that of AGNs with high X-ray luminosi ties (Stcffc n~et al.l 
l2003HUeda et al.ll2003HHasinger et all 120051. Using the 
SDSS galaxies iHeckman et a l. (2004) conclude that this 
behavior is driven by a decrease in the characteristic mass 
scale of actively accreting black holes. As the total stellar 
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Fig. 7. — R — K and the specific SFR (as derived from the stellar continua) vs. the stellar mass for both the UV-selected galaxies by 
Erb ct al. (2006a b e) and our spectroscopic X-selccted sample at 2.0 < z < 2.7. The masses f or the UV and X -selected samples are derived 
from modeling the broadband photometry and the stellar continua respectively, assuming a ISalpeterl |jl955) IMF. The galaxies in which 
AGNs are identified are indicated in black. The shaded gray areas in both diagrams are empty probably because of incompleteness effects. 
We can roughly remove the bias towards AGN with more massive black holes, by normalizing the star-formation rate by the stellar mass, 
assuming that it scales with the black hole mass. Considering a line of constant specific SFR the line emission originating from accretion 
onto the black hole should be equally easy to detect as the line emission from star formation. Thus, AGNs with the same normalized 
accretion rate should have been detected in the area between the dotted lines and to the right of the shaded region. This figure illustrates 
that black hole accretion is more effective at the high-mass end at z ~ 2.3. 



mass and black hole mass are found to correlate (e.g., 
iDesroches et al.l |2006[ ) , the stellar mass of AGN hosts 
may be expected to decrease as well, when going to lower 
rcdshift. 

To explore the behavior of the stellar mass of AGN 
hosts with redshift, we compare the AGN fractions of 
galaxies with similar normalized accretion rates at high 
and low redshift per stellar mass bin in Figure [5] The 
high-redshift galaxies are indicated by the symbols with 
error-bars. For the UV-selected galaxies we only use the 
mass-range to the right of the dashed line indicated in 
Figure [Tb, as we cannot ensure the detection of AGNs at 
lower stellar masses. The fraction peaks at M ~ 2.9 x 
10 11 M Q , but as explained in the previous subsections, 
this may be due to incompleteness effects. 

In order to compare our sample with low-redshift 
AGNs with a similar normalized accretion rate we impose 
a threshold for L[ NII ]/M«, assuming that the black hole 
mass scales with the stellar mass. For the low-redshift 
galaxies we use the same SPSS ga l axy sam ple and AGN 
classification by iKauffmann et a l. (2003b), as plotted in 
Figure [3] We required a > 3a detection for both Ha and 
[Nil]. In case [Oiii]/H/3 could not be determined, we 
selected AGNs by log([Nii]/Ha)> -0.2. We set the nor- 
malized accretion threshold to the minimum luminosity 
of all detected [N n] lines in our sample - as AGNs with 
these luminosities would have been detected - divided by 
the median stellar mass of the AGN hosts in our sample. 

At low redshift the accretion rat e is often parametrize d 
by the luminosity of [Om] (e.g.. iHeckman et al.ll2005f ). 
As [O in] is not detected for most of our galaxies, we use 
the luminosity of [Nil] instead. For AGNs in the local 
universe there is a broad positive correlation between the 
two (see Figure 3), although this correlation is different 



for LINERs and Seyfert 2s. Nevertheless, [Nil] is accept- 
able for our purpose, as we only use the line luminosity to 
determine a normalized accretion threshold to compare 
similar samples at low and high redshift, which include 
both LINERs and Seyfert 2s. 

Figure [5] shows that SDSS AGN hosts with a simi- 
lar normalized accretion rate as our high-redshift AGNs 
reside in lower mass galaxies (~ 4 x 1O 1O M0). This 
finding provides direct observational support that cos- 
mic downsizing of AGNs is related to the decrease of the 
characteristic stellar mass of their host galaxies (see also 
IHeckman et al.l 120041 ) . This conclusions is fairly robust, 
and is not strongly dependent on the threshold value 
L[mi]/M*, or whether we would use a threshold based 
on the luminosity of [O in] instead: a factor 2 or 3 differ- 
ence barely changes the stellar mass at which the AGN 
fraction peaks at low redshift. The main concern for 
using the [O in] or [N n] emission lines as indicators for 
the accretion rate is the contamination by star forma- 
tion. However, in case [Nil] at z ~ 2.3 is more contam- 
inated by star formation than at z ~ 0, our normalized 
"threshold accretion rate" would be overestimated, and 
we would select less galaxies at low redshift, and find a 
slightly lower typical stellar mass. 

Figure [8] also shows the total AGN fraction (i.e., in- 
cluding fainter AGNs) among the SDSS galaxies versus 
stellar mass. This distribution, which is dominated by 
faint AGNs, peaks at higher stellar masses. This implies 
that at similar stellar masses of the host galaxies, AGNs 
in massive galaxies at low redshift are less powerful than 
those at high redshift. These low-redshift AGNs may be 
the low luminosity descendants of the AGNs observed in 
our if-selected galaxies. 

Both these conclusions would be qualitatively similar 
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Fig. 8. — Fraction of actively accreting AGN among galax- 
ies in different mass bins at different redshifts. The red dia- 
monds present the AGN fraction for the 20 /{"-selected galaxies 
at 2.0 < z < 2.7. The purple squares present the AGN fractio n 
among the UV-selected galaxies at z ~ 2.2 by lErb et al.l {2006a). 
The rapid fall off of both histograms may be due to selection effects. 
The blue histogram is th e fraction of AGNs in the SDSS sample 
(Kauffmann ct al. 2003a) with similar normalized accretion rates 
as the AGNs in our /{"-selected sample. The green histogram are 
all AGNs in the SDSS, as defined in the text. All stellar masses 
are derived fro m stellar population modeling, and converted to a 
Salpetcr (1955) IMF. This diagram illustrates that actively accret- 
ing AGNs reside in more massive galaxies at higher redshift. 



if our AGN fraction would be slightly over- or underes- 
timated, as the galaxies in our sample span only a small 
range in stellar mass. 

4.4. AGNs and the suppression of star formation 

The decrease of the stellar ma ss of actively star- 
forming galaxies with redshift (jCowie et alj 119961 : 
Puneau et aill2005l ) may imply that the stellar mass at 
which the star formation is suppress ed also decreases over 
cosmic time. iBundv et al.l ()2006[ ) introduce a critical 
"quenching mass" , which declines significantly between 
a redshift of 1.4 and 0.5, and drops even further at low 
redshift, as the galaxy bi-modality in the SDSS breaks 
down at a stellar mass of 3 x 10 10 M Q (|Kauffmann et all 
2003a). Also, theoretical studies derive critical quench- 
ing masses, above which cooling and star form ation are 
shut down abruptly (e.g.. iCattaneo et"aill2006f h 

Due to incompleteness of our sample, especially at 
lower stellar masses, we cannot derive a critical quench- 
ing mass for the observed redshift range. Nevertheless, 
Figure[5]illustrates that in contrast with low redshift, the 
high-mass end at z ~ 2.3 consists of both star-forming 
galaxies and red, quiescent systems. Moreover, the red, 
passive systems may dominate the galaxy population be- 
yond 2.5 x 10 n M Q (see Figured]). This may imply that 
the transition from star-forming to quiescent galaxies at 
z r~j 2.3 occurs at the observed mass range. 

Remarkably, at both low and high redshift, the fraction 
of actively accreting AGNs may peak at a similar stellar 
mass range at which the suppression of the star formation 



seems to occur (see also lHeckman et a l. 2004). This may 
imply that the critical quenching mass tracks the evolu- 
tion of the typical stellar or black hole mass at which ac- 
tive accretion takes place. This is illustrated in Figure [9l 
which shows the rest- frame u — g color and D n (4000) ver- 
sus the stellar mass for both our AT-selected sample and 
the 28 000 low-redshift SDSS galaxies with estimated co- 
moving distan c es in the range 10 < d < 150 Mpc h _1 by 
iBlanton et~aT1 ()2005f ) (this is a sub-sample of the SDSS 
galaxies presented in Figures [3] and [5J) . D n (4000) mea- 
sures the ratio of the average flux density F v in the 
bands 4000-4100 A and 3850-3950 A around the break 
(Balo gh~et al.l li999). The actively accreting low-redshift 
AGNs are identified by the orange contours. 

The AGN hosts at low and high redshift seem to have 
similar D„(4000) and rest-frame u — g colors. However, 
as explained in Section |3~31 AGNs in star forming galax- 
ies may have been missed, and this would lower the 
typical £>„(4000) and rest-frame u — g colors of AGN 
host galaxies. We note however, that the AGN classi- 
ficat ion in the SPSS prob ably suffers from a simila r ef- 
fect. ISanchez et ail (j2004D and lNandra et alj (|2006l) find 
a similar location for AGNs in color space - in the "green 
valley" between the red sequence and blue cloud - for op- 
tically and X-ray selected AGNs at 0.5 < z < 1.1 and 
0.6 < z < 1.4 respectively. Although the AGN host 
galaxies at both low and high redshift and the evolved 
galaxies at high redshift have colors comparable to the 
red-sequence galaxies at low redshift, their 4000 A breaks 
are much smaller (Fig. [9]). The red colors of these galax- 
ies mainly reflect their strong Balmcr breaks (instead of 
the 4000 A break), indicative of a post-starburst SED. 

In summary, our results may suggest that the suppres- 
sion of star formation is related with an AGN phase. 
Different models based on the CDM theory use vari- 
ous types of AGN fee dback to stop t he sta r formation 
in massive galaxies. iHopkins et al.l (|2006h propose a 
quas ar mode which is v ery effective at high redshift, 
while ICroton et ail (|2006l ) introduce a radio mode which 
is more efficient at low redshift and can also account for 
the maintenance of "dead" galaxies. Our result might 
be more cons i stent with the mechanism proposed by 
I Hopkins et alJ (|2006[ ). as the radio mode b ylCroton et alJ 
( 2006) is not very effective at z > 2. In the Ho pkins et al.l 
(2006) picture, the observed emission lines may be the 
fading activity of a preceding quasar phase. 

5. SUMMARY AND CONCLUSIONS 

In order to understand the formation of massive galax- 
ies at high- redshift, we are conducting a NIR spectro- 
scopic survey of massive galaxies at z ~ 2.5. Our K- 
selected galaxy sample contains 20 spectroscopically con- 
firmed galaxies at 2.0 < z < 2.7. For nine of these galax- 
ies no emission lines are detected, indicating that the star 
formation in these s ystems is already strongly suppressed 
(|Kriek et a l. 2006b). In this paper we focus on the emis- 
sion line galaxies in the sample, and attempt to constrain 
the main origin of the ionized emission. Remarkably, we 
find that at least four of the 20 galaxies in our sample 
may host an AGN. These four AGN candidates are iden- 
tified by a combination of indicators, such as the emis- 
sion line ratios, compactness of the line emission, and 
the star formation rate of galaxies as derived from the 
stellar continua. Furthermore, there are three additional 
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Fig. 9. — Rest-frame u — g color and D n (4000) vs. stellar mass for the i\-selected 2.0 < z < 2.7 galaxy sample. Both properties (u — g and 
D n (4000)) are derived from the best fits to the rest-frame optical continuum spectra in combination with the rest-frame UV broad-band 
photometry. The errors are derived from the 68 % best-fits of 2 00 M onte Carlo simulations. T he gray scale presents th e low-redshift 
catalogue of galaxies extracted from the SDSS by Blanton ct al. (2005). The photometric masses (Kauffmann ct al. 2003a) are corrected 
for the difference in assumed IMF with the if-selected sample. The actively accreting AGNs within the SDSS sample, as defined in § 14.31 
are represented by the orange contours. The evolution of a stellar population for two different SFHs with different decaying times (t) are 
drawn, assuming a formation redshift of 2.8. The ages are indicated in Gyr. These diagrams may suggest that both at low and high redshift 
the AGN activity peaks during the post-starburst phase, when a galaxy is transforming from a star-forming into a quiescent galaxy. 



galaxies with [N n] /Ha ratios that may be indicative for 
AGNs. However, this hypothesis is not supported by the 
other diagnostics (compactness of line emission, Wn a , 
MIR spectral shape), and other ionization mechanisms, 
like starburst driven winds also seem plausible for these 
galaxies. 

This work has several direct implications for high red- 
shift studies. First, as for a significant part of the mas- 
sive galaxies at 2.0 < z < 2.7 the line emission is not 
dominated by pure photoionization processes from H II 
regions, the Ha luminosity is not a good diagnostic for 
the current star formation rate. Second, none of the 
AGN candidates has a counterpart in the 100 ks (for 
SDSS1030) or 250 ms (for ECDFS) X-ray images. This 
implies that optical emission line ratios provide a com- 
plementary approa ch to identify (relativ ely low luminos- 
ity) AGNs (see also Heckman et al. 2005, for low-redshift 
AGNs). 

The stellar populations of the four AGN hosts range 
from evolved to star forming, with a median absolute and 
specific SFR of 9M© yr _1 and 0.04 Gyr -1 respectively. 
However, this result may be biased, as AGNs are easier 
to identify in evolved galaxies. 

Combining our if-selected sample with a UV-selected 
sample in the same redshift range, that spans a much 
broader range in stellar mass, we find that black- hole 
accretion is more effective at the high-mass end of the 
galaxy distribution (-2.9 x 10 n M Q )atz~2.3. This 
result may be partly due to selection effects, as both 
samples suffer from incompleteness, and AGNs with the 
same normalized accretion rate are easier to identify in 
more massive systems. 

The AGNs in low redshift SDSS galaxies with the same 
stellar mass as our if-selected galaxies are less luminous 
(as implied from their [N n] line luminosities) than their 



high-redshift analogs. These AGNs may be the low- 
luminosity descendants of those we detect at z > 2. How- 
ever, AGNs with similar normalized (to A/* — Mbh ) 
accretion rates as the z > 2 AGNs reside in less massive 
galaxies (~ 4 x 1O 1O M ) at low redshift. This is di- 
rect observational evidence for downsizing of AGN host 
galaxies. It will be interesting to see when the AGNs in 
the z ~ 2.3 massiv e galaxies switch off. We note that 
iBarger et ail (|2005l) show that AGNs in more optically 
luminous, and thus probably more massive host galaxies 
are switching off between z=l and the present time. 

In contrast to what we see at low redshift, the mas- 
sive galaxy population at z — 2.3 is very diverse. While 
the red sequence is starting to build up, half of the mas- 
sive galaxies are still vigorously forming stars. The wide 
spread of properties observed in massive galaxies at this 
redshift range, suggests that we sample different evolu- 
tionary stages of the processes that transform a blue 
galaxy into a red one. Thus, the transition may take 
place at the same mass range at which the AGN fraction 
may peak. This is similar to what we see at low redshift: 
the fraction of AGNs with the same normalized accre- 
tion rate as those observed at high redshift peaks at the 
same stellar mass range at which the color bi-modality 
breaks down. Furthermore, the rest-frame u — g colors 
and Z?„(4000) of the AGN host galaxies at low and high 
redshift span the same range, and show that actively ac- 
creting AGNs mainly reside in post-starburst galaxies. 
This may suggest that the suppression of the star for- 
mation is related to an AGN phase, and that the typical 
stellar mass scale at which this suppression occurs de- 
creases with redshift. Our results are qualitatively con- 
sistent with the AGN feedback m echanism as propose d 
by several theoretical studies (e.g.. iHopkins et aT1 r2006). 
However, the detailed mechanism is not well understood, 
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and perhaps our observations can help to constrain the 
theory. 

The fact that three of the AGN host galaxies have a 
very low SFR may suggest that feedback processes are 
important. However, we do not have direct evidence that 
the suppression of the star formation is caused by the 
AGNs, as the co-evolution of star formation and AGN 
activity can simply be explained by the fact that both 
the AGN and star-formation are fueled by a large gas 
supply. In this picture, the galaxy simply runs out of gas, 
leading to a fading of the AGN and the star formation at 
approximately the same time. To clarify this situation it 
is necessary to determine whether massive star-forming 
galaxies in the studied redshift range host AGNs as well, 
and how their accretion rates relate to the accretion rates 
of the AGNs in the post-starburst galaxies. 

We also note that the AGN classification is somewhat 
uncertain for each individual object, and we may 
have under- or over-estimated the AGN fraction. Our 
results are qualitatively similar if one or two galaxies 
were misclassified (in either direction). However, our 
conclusions obviously do not hold if none of the galaxies 
in the sample host an AGN. As the physical origin of 
LINER emission is not well understood even at low 



redshift, we cannot exclude the possibility that the 
observed emission line ratios at high redshift originate 
from physical processes other than black hole accre- 
tion. Deeper data than presented in this paper for a 
significantly larger sample, in combination with optical 
spectroscopy, deep X-ray data and MIR photometry 
should address this concern and shed new light on the 
question what processes are responsible for transforming 
a star-forming galaxy into a red, passive system. 



We thank the referee for thorough and helpful com- 
ments, and Shanil Virani and Carie Cardamone for 
their help in obtaining the X-ray upper limits. This 
research was supported by grants from the Nether- 
lands Foundation for Research (NWO), and the Leids 
Kerkhoven-Bosscha Fonds. Support from National Sci- 
ence Foundation grant NSF CAREER AST-0449678 is 
gratefully acknowledged. DM is supported by NASA 
LTSA NNG04GE12G. EG is supported by no. AST- 
0201667, an NSF Astronomy and Astrophysics Postdoc- 
toral Fellowship ( AAPF) . PL acknowledges support from 
Fondecyt Grant no. 1040719. 



REFERENCES 



Alonso-Herrero, A., et al. 2006, ApJ, 640, 167 

Baldwin, J., Philips, M., & Terlevich, R. 1981, PASP, 93, 5 

Balogh, M.L., Morris, S.L., Yee, H.K.C., Carlberg, R.G., & 

Ellingson, E. 1999, ApJ, 527, 54 
Barger, A. J., Cowie, L.L., Mushotzky, R.F., Yang, Y., Wang, W.- 

H., Steffen, A.T., & Capak, P. 2005, AJ, 129, 578 
Bonnet, H. et al. 2004, The ESO Messenger 117, 17 
Bower, R.G., et al. 2006, MNRAS, 370, 645 
Blanton, M.R., et al. 2005, AJ, 129, 2562 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Bundy, K, et al. 2006, ApJ, 651, 120 
Calzetti, D. 1997, AJ, 113, 162 

Calzetti, D., Armus, L., Bohlin, R.C., Kinney, A.L., Koornheef, J., 

& Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Cattaneo, A., Dekel, A., Devriendt, J., Guiderdoni, B., &; Blaizot, 

J. 2006, MNRAS, 370, 1651 
Cecil, G., et al. 2000, ApJ, 536, 675 

Cid Fernandes, R., Mateus, A., Sodre, L., Stasihska, G., & Gomes, 

J.M., 2005, MNRAS, 358, 363 
Cowie, L.L., Songaila, A., Hu, E., & Cohen, J.G. 1996, AJ, 112, 

839 

Croton, D.J., et al. 2006, MNRAS, 365, 11 
Daddi, E., et al. 2005, ApJ, 600, L127 
Daddi, E., et al. 2004, ApJ, 617, 746 

Desroches, L.-B., Qua taert, E., Ma, C.-P ., & West, A. A. 2006, 

MNRAS, submitted l |astro-ph/0608474} 
Dopita, M.A., & Sutherland, R.S. 1995, ApJ, 455, 468 
Eisenhauer, F. et al. 2003, SPIE 4841, 1548 
Elias, J.H., et al. 2006, SPIE 6269, 139 

Erb, D.K., Shapley, A.E., Pettini, M., Steidel, C.C., Reddy, N.A., 

& Adelberger, K.L. 2006, ApJ, 644, 813 
Erb, D.K., Steidel, C.C., Shapley, A.E., Pettini, M., Reddy, N.A., 

k, Adelberger, K.L. 2006, ApJ, 646, 107 
Erb, D.K., Steidel, C.C., Shapley, A.E., Pettini, M., Reddy, N.A., 

k, Adelberger, K.L. 2006, ApJ 647, 128 
Fabian, A.C., et al. 2003, MNRAS, 344, L43 

Fabian, A.C., Sanders, J.S., Taylor, G.B., Allen, S.W., Crawford, 
C.S., Johnstone, R.M., & Iwasawa, K. 2006, MNRAS, 366, 41 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Filho, M.E., Fraternali, F., Markoff, S., Nagar, N.M., Barthel, P.D., 
Ho, L.C., & Yuan, F. 2004, A&A, 418, 429 

Forster Schreiber, N.M., et al. 2004, ApJ, 616, 40 

Forster Schreiber, N.M., et al. 2006a, AJ, 131, 1891 

Forster Schreiber, N.M., et al. 2006b, ApJ, 645, 1062 

Franx, M. et al. 2003, ApJ, 587, L79 

Gawiser, E. et al. 2006, ApJS, 162, 1 



Gebhardt, K., et al. 2000, ApJ, 539, L13 
Genzel, R., et al. 2006, Nature, 442, 786 
Giacconi, R., et al. 2002, ApJS, 139, 369 

Granato, G.L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 

2004, ApJ, 600, 580 

Hasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Heckman, T.M., Kauffmann, G., Brinchmann, J., Chariot, S., 

Tremonti, C, & White, S.D.M 2004, ApJ, 613, 109 
Heckman, T.M., Ptak, A., Hornschemeier, A., & Kauffmann, G., 

2005, ApJ, 619, 35 

Hopkins, P.F., Hernquist, L., Cox,T.J., Robertson, B., & Springel, 

V. 2006, ApJS, 163, 50 
Juneau, S., et al. 2005, ApJ, 619, L135 

K ang, X., Jing, Y .P., & Silk, J. 2006, ApJ in press 

jastro-ph/ 0601685 ) 
Kauffmann G., et al, 2003, MNRAS, 341, 33 
Kauffmann G., et al. 2003, MNRAS, 346, 1055 
Kennicutt, R.C. 1998, ARA&A, 36, 189 

Kewley, L.J., Dopita, M., Sutherland, R., Heisler, C, & Trevena, 

J. 2001, ApJ, 556, 121 
Kewley, L.J., Groves, B., Kauffmann, G., & Heckman, T. 2006, 

MNRAS, 372, 961 
Kriek, M., et al. 2006a, ApJ, 645, 44 
Kriek, M., et al. 2006b, ApJ, 649, L71 
Labbe, I., et al. 2005, ApJ, 624, L81 
Li'pari, S., et al. 2004, MNRAS, 355, 641 
Nandra, K, et al. 2006, ApJ, in press l|astro-ph /0607270 ) 
Nesvadba, N.P.H., Lehnert, M.D., Eisenhauer, F., Gilbe 

Tecza, M., & Abuter, R. 2006 ApJ, 650, 693 
Papovich, C, et al. 2006, ApJ, 640, 92 
Quadri, R., et al. 2007, ApJ, 654, 138 

Reddy, N.A., Erb, D.K., Steidel, C.C., Shapley, A.E., Adelberger, 

K.L., & Pettini, M. 2005, ApJ, 633, 748 
Reddy, N.A., et al. 2006, 644, 792 

Rubin, K.H.R., van Dokkum, P.G., Coppi, P., Johnson, O., Forster 

Schreiber, N.M., Franx, M., & van der Werf, P. 2005, 613, L5 
Salpeter, E.E. 1955, ApJ, 121, 161 
Sanchez, S.F. et al. 2004, ApJ, 614, 586 
Schawinski, K, et al. 2006, Nature, 442, 888 

Shapley, A.E., Coil, A.L., Ma, C.-P., & Bundy, K. 2005, ApJ, 635, 
1006 

Steffen, A.T., Barger, A.J., Cowie, L.L., Mushotsky, R.F., & Yang, 

Y. 2003, ApJ, 596, L23 
Steidel, C.C., Hunt, M.P., Shapley, A.E., Adelberger, K.L., Pettini, 

M., Dickinson, M., & Giavalisco, M. 2002, ApJ, 576, 653 



Gilbert, A. 



16 



Kriek et al. 



Steidel, C.C., Shapley, A.E., Pettini, M., Adelberger, K.L., Erb, 

D.K., Reddy, N.A., Hunt, M.P. 2004, ApJ, 604, 534 
Stern, D., et al. 2006, ApJ, submitted (astro-ph/0608603 ) 
Tremonti, C.A. et al. 2004, ApJ, 613, 898 

Ueda, Y., Masayuki, A., Ohta, K., & Miyaji, T. 2003, ApJ, 598, 
886 

van Dokkum, P.G. 2001, PASP, 113, 1420 
van Dokkum, P.G. et al. 2004, ApJ, 611, 703 

van Dokkum, P.G., Kriek, M., Rodgers, B., Franx, M., & Puxley, 
P. 2005, ApJ, 622, L13 



van Dokkum, P.G., et al. 2006, ApJ, 638, L59 
van Dokkum, P.G., & van der Marel, R.P., 2006, ApJ, 655, 30 
Veilleux, S., & Osterbrock, D. 1987, ApJS, 63, 295 
Virani, S.N., Treister, E., Urry, CM., & Gawiser, E. 2006, AJ, 131, 
2373 

Weedman, D., et al. 2006, ApJ, 653, 101 
Wuyts, S., et al. 2006, ApJ, 655, 51 



